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Abstract 

The nuclear enhancement observed in inelastic photoproduction of J/ ^ should not 
be interpreted as evidence for an increased gluon density in nuclei. The nuclear sup- 
pression of the production rate due to initial and final state interactions is calculated 
and a novel two-step color exchange process is proposed, which is able to explain the 
data. 
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1. Introduction 

Measuring inelastic photoproduction of J/\I/ off heavy nuclei, J/\E'X, the EMC 

collaboration |^ and later the NMC collaboration have observed an enhancement of the 
production rate compared to a proton or a light nuclear target: the experimental ratio of 
the cross sections per nucleon on tin to carbon exceeds the value one by about 10%. If this 
enhancement is fully attributed to an enhancement of the gluon density - as the authors of 
experiment do - the gluon density in heavy nuclei should be larger by a similar value of 
a few percent, which is perfectly compatible with the theoretical expectations. 

This conclusion is dangerous, however, since it neglects initial and final state interactions 
due to the presence of the other nucleons. The data from the E772 experiment |§ for J/^ 
production in proton-nucleus collisions show a suppression, which may be parameterized by 
A"*^, where e = 0.08. There is no reason, why this kind of initial/final state interaction should 
not be present in photoproduction. Therefore, the observed enhancement in the Sn/C ratio 
in 7 J/\l/ production of the order of 10% has to be seen in the light of an expected 
20% suppression due to initial/final state interaction. It is theoretically unacceptable to 
claim that the discrepancy of 30% is attributed to a modified gluon distribution. Other 
mechanisms have to be looked for. This is the aim of the present paper, which is organized 
as follows. 

In the next section the triple-Regge approach to inelastic photoproduction of charmo- 
nium is developed. Since the charmonium is produced diffractively via Pomeron exchange, 
we conclude that at Xi ^ 1 the reaction is sensitive to the double-gluon distribution function 
(similarly as in elastic photoproduction), rather than to the single gluon density Then, 
in section 3, the suppression factor for inelastic photoproduction of J/^ off nuclei is calcu- 
lated using Glauber's eikonal approximation for the final state interaction. The formulas 
are derived for the first time. Since Glauber's approximation assumes lack of formation time 
we consider the color-singlet model for J/\E' photoproduction in section 4. We demonstrate 
that contrary to the wide spread opinion the time interval for evaporation of color is quite 
short, about 1 fm in the energy range of the NMC experiment. Therefore Glauber's model 
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can be used in this case as well. Inelastic shadowing corrections to Glauber's approximation 
are considered in section 5. We demonstrate that this corrections make the nuclear medium 
substantially more transparent for charmonium: the effective absorption cross section turns 
out to be only about 60 % of that for J/^*. In section 6 we present the novel mechanism for 
diffractive photoproduction of charmonium off nuclei: the photon interacts with a bound 
nucleon via color-exchange and produces a cc pair in a color-octet state. The pair prop- 
agates through nuclear matter and looses energy via hadronization. Then, with another 
nucleon the cc pair exchanges color and converts to a color-singlet state. Such a process 
leads to J/^^'s with xi < 1. Since this mechanism relies upon two interactions at least, its 
cross section may have an 74"-dependence with a > 1. Interplay of this mechanism with 
the standard multiple scattering one nicely explains the NMC data. 

2. Glauber theory for inelastic J/^ photoproduction off nuclei 

The inelastic photoproduction cross section of a J/^' on a nucleus can be represented in 
the form 

dxi * dxi 

In this section we derive the nuclear suppression factor S^"^ within Glauber theory. There 
are two types of amplitudes which contribute to the inelastic photoproduction of a J/^' off 
nuclei: 

I) The incoming photon does not interact in the nucleus until the point with coordinates 
{b,z) (the impact parameter and the longitudinal coordinate, respectively), where it inter- 
acts inelastically, — > J/^X, and produces the J/^ with the detected final momentum, 
which leaves the nucleus without inelastic interactions; 

II) prior the inelastic interaction at the point (6, z) the photon produces elastically a 
J/^ at the point (&, -^i), — > J/^N, with zi < z. The J/^ propagates through the 
nucleus without interaction up to the point {b,z), where it produces the final J/^ in the 
reaction J/^'iV J/vfX. 
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In both cases the longitudinal momentum transfer at the point (6, z) is so large, qin ~ 
(1 — Xi)z/, that only those amplitudes interfere, which correspond to the inelastic production 
of the J/\l/ with the desired momentum on the same nucleon. 

On the other hand, the longitudinal momentum transfer in elastic photoproduction at 
the point (6, Zi) may be small at high energies, 

del = . (2) 

Therefore the amplitudes corresponding to different coordinates Zi may interfere. 

According to our classification of the amplitudes I and II the suppression factor S^'^ has 
three terms, S'^'^ = (5*^^)1 + (5^^)12 + {Sl'^)2, where the first and the third terms correspond 
to squares of the amplitudes I and II, respectively, and the second one is the interference 
term. 

The amplitude I squared provides the following suppression factor 

00 

{Sr), = jjd%ldz pib,z) exp [-af/ T,(6)] . (3) 
—00 

Here Tz{b) = dz'p{b, z'), where p{b, z) is the nuclear density normalized to A. The nuclear 
thickness function is T{b) = Tz^^^oib)- 
The interference term has a form 
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This term has a negative sign because the amplitude II plays the role of an absorptive 
correction to amplitude I. 

To calculate the contribution of the amplitude II squared we can use the results of ref. 
|[T0| , PI for quasielastic photoproduction of vector mesons on nuclei, which we should apply 
to production of the J/\l/ prior the point {b,z). However, as different from |n[], we 



should include both coherent (small qt) and incoherent (large qt) interactions of the J/\E' 
in the nucleus. The result reads, 
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At low energies, when q^i 3> 1/Ra only (^Sq/^^^ and the first term in curly brackets in 
survive and lead to the nuclear suppression factor 
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At high energies, in the limit Qei — > the sum of expressions (§)-(§) simplifies very much, 

SlX.<<./n. = jld%Tib)e-'^l''-^'K (7) 

Using the smallness of af^^ ~ 5 mb we can expand (D as S^'^lg^^^i/R^ 1 — |cr*f^ (T), 
where the mean nuclear thickness (T) = J d% T'^{b). In the approximation eq. (^ we 
have kept only the first order of the expansion parameter o"*f^(T), and neglected cr*^(T), 
which is of the second order. To the same order the low-energy limit eq. (P) leads to 

Thus, the nuclear suppression factor at high energy turns out to be smaller than that at 
low energy. This result has a natural space-time interpretation: at high energies a photon 
may interact strongly via its hadronic fluctuations. The lifetime of the cc fluctuation of 
the photon, called coherence time, is tc ~ l/^ez- A fluctuation created long in advance 
propagates and attenuates through the whole nuclear thickness, as is explicitly exposed 
in (1^. On the other hand, at low energy the photon converts to a cc inside the nucleus 
just prior to the J/\E' production. So the cc pair propagates only through about a half 



of the nuclear thickness. Such a contribution corresponds to {S^'^)i given by (j^). There 
is, however, a possibihty to produce J/\l/ "elastically" in —>■ J/'^N prior the inelastic 
rescattering Jj^N Jj^X. This correction is suppressed by factor crfi^ /afj^ and is also 
included in eq. (|^). 

Expressions (H) - give the correct extrapolation between the low-energy eq. (j^) and 
the high-energy eq. limits and are quite complicated. If one uses the same approximation 
afj^ (T) <^ 1 and expands eqs. (0) - (H) in this small parameter up to first order (compare 
with quasielastic photoproduction of J/\E' off nuclei [12, 13||), one finds 
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where the effective longitudinal nuclear formfactor squared is defined as 
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One can approximate Fj^lqei) quite well using a Gaussian parameterization instead of the 
more realistic Woods-Saxon form for the nuclear density and obtains F^gez) = exp(— i?^ Q^i)/^-, 
where R\ is the mean squared charge radius of the nucleus Note that S'^^ is indepen- 
dent of Xi, therefore this mechanism of nuclear suppression keeps the same Xi-distribution 
of the cross section as on a free nucleon and cannot account for the data (cf. section 4). 

There is another possibility within the standard Glauber approach to produce the J/\l/ 
in two-step process on a nucleus: the photon produces the J/\l/ inelastically on one nucleon 
with a subsequent diffractive scattering of the charmonium on another nucleon. Inelasticities 
of these two interactions should be adjusted to provide the desired xi for the final J/\l/. The 
contribution of this mechanism to the cross section is evaluated in the Appendix and a value 
of less than 2% is found, which we neglect compared to the main mechanism considered 
above. 



3. Triple- Regge formalism for J/^ photoproduction. 
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As soon as the nuclear suppression factor Sl'^ is calculated within Glauber approxima- 
tion, one can predict the nuclear cross section provided that the J/\l/ photoproduction cross 
section on a free nucleon is known. We review briefly in this and the next sections the well 
accepted mechanisms in order to justify the usage of Glauber model and to get a reasonable 
parameterization for the cross section. 

In the standard triple- Regge formalism (see for instance 0) the diffractive inelastic 
photoproduction of J/\E' is described by the triple- Pomeron graph, depicted in Fig. 1. 




y 



Figure 1: The triple- Pomeron graph for the cross section of the reaction 
7A^ J/'^X. The dashed line shows that only the absorptive part of 
the amplitude is included. 



The corresponding differential cross section reads 

ia(,N -,J/^X) ^ ^ IT'^Sf, . (10) 

dxi dt ^ ' ' ^ (1 - xi)2"p(*)-i ^ ^ 

Here Xi and t are the relative share of the photon light-cone momentum carried by the J/\& 
and the 4-momentum squared of the produced J/\l/, respectively. ap{t) is the Pomeron 
trajectory and Bfj' is the slope parameter for inelastic J — p collisions. The factor 
G{'~fN —>■ J/'^X) includes the triple- Pomeron coupling as well as the P'-j^! and the PNN 
vertices. It can be evaluated using the data from hadronic interactions, assuming factoriza- 
tion of the Pomeron. 
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G{^N J/^X) = G{pp ^ pX) ""^^P ^J^^^^ . (11) 

The momentum transfer squared t is rather large in J/\E'X, = (M| + 

XiQ'^){l — xi)/xi, where Xi = p^/p^ (~ xp at Xi 1) is the fraction of the photon hght- 
cone momentum carried by the J/"^. \tmin\ may reach a few GeV^ dependent on xi. The 
Pomeron trajectory ap{t) is known to approach the value of 1 at high \t\ 0], therefore we 
take ap{t) = 1. 

Combining all the factors we get an estimate for the Xi-distribution of the inclusive 
photoproduction of J/^' on a proton, integrated over transverse momenta 

da(^N ^ J/M/X) ^(7P ^ J/'^P) Gjpp ^ pX) 

dx, BlfBl^ <f l-x, ■ ^ ^ 

The effective triple- Pomeron constant G{pp pX) = 3.2 mb/GeV is fixed by an analysis 
of data on pp — > pX [^]. The elastic pp cross section and the slope parameter are taken 
as cTgf = 8 mb and B^f ^ 10 GeV^"^. The slope of the pr-distribution of elastic J/\l/ 
photoproduction is assumed Bfi^ ^ B^f /2, since the slopes are related to the nucleon and 
the J/\l/ form factors respectively. On the other hand, in the inelastic photoproduction 
of J/\E' the nucleon is destroyed and its formfactor does not contribute to the effective 
slope (the same as in the case of pp — » pX reaction, where the slope Bf^ ^ 4 GeV^^ 
Therefore BfJ^ is only due to the charmonium formfactor. It was measured in the EMC 
experiment § to be bZ^ = 0.58 ± 0.07 GeV'^, but the NMC experiment § led to Bfj' = 
0.23 ±0.02 GeV~^. For a mean photon energy (u) ^ 100 GeV, the elastic photoproduction 
cross section is a{'~fp J/'^p) ^ 12 nb according to the EMC measurements []^, but is 
^ 27 nb according to the NMC results In view of the uncertainties in the measured 
values of the elastic photoproduction cross section and the inelastic slope we fix them at a 
middle value a^-fp J/^p) = 20 nb and Bfj' = 0.45 GeV-"^. 

The the cross section formula (0) is calculated with the parameters given above and is 
compared with the data from the NMC experiment in Fig. 2. A good agreement is found 
both in shape and normalization. 
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Figure 2: The data from the NMC experiment P| for inelastic J/\l/ pho- 
toproduction on protons. The sohd curve shows predictions based on the 
triple-Pomeron phenomenology. 

Note that the triple-Pomeron is known to dominate in the inclusive cross section of pp 
interaction only at Xi > 0.9 0. The nondiffractive RRP triple- Reggeon graphs are impor- 
tant at smaller Xi. OZI rule suppresses the RRP term for J/\E' interactions, therefore the 
triple Pomeron contribution may be important down to lower Xi. Nevertheless, the triple- 
Pomeron phenomenology does not necessarily describe the data for xi < 0.8. Particularly, 
the longitudinal formfactor of the 7 — J/^ vertex may substantially diminish the cross 
section at mid Xi. Yet we need a parameterization of the data for further applications, and 
Fig. 2 shows that eq. (|T2|) reproduces the data well down to Xi = 0.6. 

Usually the inelastic photoproduction of J/\l' is treated perturbatively, within the so 
called color-singlet model 0]. In that approach the cross section is predicted to be pro- 
portional to the low-x gluon density distribution in the proton. Although the Pomeron by 
definition includes all possible kinds of perturbative graphs, but the ones of the color-singlet 
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model are not the leading contributions at Xi ^ 1. Therefore the spectrum of J/\l/ predicted 
by the color singlet model does not peak at a;i ^ 1, while the triple- Pomeron contribution 
does. The dominant mechanism comes from the ladder-type graphs with excitation of the 
proton. Correspondingly, the cross section is proportional to the double-gluon density, which 
is approximately equal to the single-gluon density squared. 

We conclude that for xi 1 the inelastic photoproduction of J/\E' cannot be used as a 
probe for the single-gluon distribution, but it might be possible at mid Xi-values. 

4. Color-singlet model 

A popular approach to inelastic diffractive photoproduction of J/\l/ on a nucleon target 
is the color-singlet model (CSM) It is somewhat similar to the mechanism suggested 
in section 6, except that the color-octet cc pair becomes colorless via radiation of a gluon, 
rather than by means of a color-exchange rescattering. 

Radiation of a gluon takes some time, which may be important if the process takes 
place in a nuclear environment. Indeed, rescatterings of the color-octet cc pair may induce 
additional energy loss, and may result in an effective A-dependence. Therefore we need to 
understand the space-time pattern of gluon radiation in the CSM. 

Sometimes the production of a charmonium in CSM is associated with a two-step space- 
time development: first the color-octet cc pair is produced, then it propagates and finally 
converts into a color-singlet state via radiation of a gluon. The second stage is usually 
assumed to take a long time. 

We do not think that such a treatment is possible in view of a principal quantum- 
mechanical uncertainty: analogous to the standard electromagnetic bremsstrahlung one 
cannot say whether the gluon/photon has been radiated before or after the interaction. 
Both Feynman graphs contribute to the cross section (actually both are taken into account 
in perturbative calculations [Q). One can understand this uncertainty as related to the 
lifetime of a fiuctuation on the initial particle containing the radiated gluon and the color- 
octet cc pair. This lifetime given by the energy- denominator equals to 
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where M is the effective mass of the fluctuation, 



M2 = ^ + . (14) 

Xi Xi[l — Xi) 

Taking the mass of the cc pair M^c ~ Mx^ we get the following expression for the radiation 
time 

_ 2z/xi(l-xi) 

fc2 +M|(l-a;i) • ^'^^ 

This radiation time is only about 1 fm in the kinematical region of the NMC experiment. 

One can safely neglect an influence of nuclear medium on the process of hadronization for 

such a short time interval. Therefore the multiple scattering approach developed in sections 

2 can be applied in the case of CSM. 

Note that even at very high energies, when becomes long, the effects of induced 

radiation by the color-octet cc pair during radiation time can be neglected since 2Kj„rf/M| <^ 

1, where Kind < 1 GeV^ is the density of energy loss via induced radiation. 

5. Beyond Glauber approximation 

Considering photoproduction of J/'^ in the quark representation one understands that 
a colorless cc wave packet, rather than a J/'^, is produced and propagates through the 
nucleus, but the J/'^ may be formed far away from the nucleus if the energy is high. Due to 
quark-hadron duality one can consider the same effects taking into account photoproduction 
of higher diffractive excitations of the charmonium, \E'', ^" etc., as well as all diagonal and 
off-diagonal diffractive rescatterings of these states in nuclear matter. The corresponding 
corrections are know as inelastic shadowing [T^. The simplest, but quite accurate two- 
coupled-channel approximation including J/\l' and \1'' was considered in |T^. It was found 
that the the cc wave packet produced in pp interaction is such a combination of J/'^ and \E'', 
which is quite close to the eigenstate of interaction having the minimal possible interaction 
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cross section. This fact explains particularly why J/\l/ and \1/ experience similar nuclear 
attenuation (cf. the E772 experiment 

We assume that the ratio of the photoproduction amplitudes 



where / is the operator of diffractive amplitude, is close to one. In pp scattering the 
corresponding value of Rpp is measured to be 0.5 (see in This assumption on R^p can 

be motivated by the similarity of the cc components of the wave functions of a photon and 
of a gluon. 

The solution of evolution equation for the cc wave packet propagated through nuclear 
matter, which incorporates with both effects of the coherence and formation lengths will be 
presented elsewhere [0, but some numerical results are published in ||T^. Here we use the 



approximation of small af^f (T) <^ 1. Then one can modify the Glauber model expression 



as follows ||T7| 



Sl^^^l- (JeffTeff ~ exp{-aeffTeff) , (17) 

where 



Te// = i(T)[l + Fi(gJ] (18) 
is according to eq. (W) the effective nuclear thickness covered by the charmonium, which is 



controlled by the coherence length Ic ~ l/o'ei in|' ^^'^ 



^eff = cTl[l + eRFl{qf)] (19) 

is the effective absorptive cross section for charmonium, which is controlled by the formation 
length for the charmonium wave function. If ^ 1/^'/, where 

_ M|,-M| 
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is the longitudinal momentum transfer in the off-diagonal reaction "^N ^ 
The new parameter e in eq. ( |l9l) is 



\ 1^1 / (^21) 



mm 



and is estimated in |16] in the oscillator approximation, e = —y2/3. Remarkably, with 
these parameters at high energies cTg// ~ 0.6 crfj^, i.e. the J/"^ attenuates in nuclear 
matter substantially less than predicts Glauber model. This may be understood as a natural 
consequence of color transparency. 

We calculate ratio of production rates of J/\I' for tin to carbon using eq. (|1^. The 
results for u = 70 and 100 GeV are plotted in Fig. 3 in comparison with the data from the 
NMC experiment [H. We see that the discrepancy with the inelastic photopro duct ion data 
(xi < 0.85) is so large that it is difficult to explain it by a few percent enhancement of the 
gluon density in heavy nuclei. One has to look for another explanation of the EMC-NMC 
effect. 

6. Novel mechanism of J/^ photoproduction in nuclear matter 

According to the arguments in section 3 the photoproduction of J/\& on a nucleon pro- 
ceeds via Pomeron exchange, which in the language of perturbative QCD can be interpreted 
as a colorless gluonic exchange, containing at least two gluons in the t-channel. In a nuclear 
interaction this colorless exchange can be split into two color-exchange interactions with 
different nucleons, keeping the beam and the nuclear remnants in the final state colorless 
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The same physics can be represented as a consequence of unitarity: the imaginary part 
of the elastic forward scattering amplitude of J/\l/ on a nucleon is a product of an inelastic 
J N amplitude with itself time-conjugated (up to small elastic and diffractive corrections). 
In the case of a nuclear target one has to consider also products of two inelastic amplitudes 
on different nucleons. 

Thus, the J/\l/ photoproduction on a nucleus can result from the photoproduction on a 
bound nucleon of a cc pair in a color-octet state, which propagates through nuclear matter 
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Figure 3: The data from the NMC experiment 0] for the ratio of inelas- 
tic J/^ photoproduction rates on tin to carbon, measured with 200 GeV 
(squares) and 280 GeV (circles) muons. The thin solid and dashed curves 
correspond to the Glauber-Gribov mechanism for nuclear suppression 
of J/\I' and V = 100 and 70 GeV, respectively. The thick solid and 
dashed curves show predictions based on the proposed novel mechanism 
of diffractive interaction with a nucleus aX v = 100 and 70 GeV , respec- 
tively. For the origin of the dividing line at xi = 0.85 see text. 



loosing energy via hadronization until another color-exchange interaction inside the nucleus 
converts the pair back to a colorless state. Actually, this is a specific QCD mechanism of 
diffractive excitation of the nucleus. The energy of the J/\E' (and the excitation energy of 
the nucleus) is less than the photon energy by the amount lost via hadronization by cc pair 
while it is in the color octet state. We assume a constant energy loss per unit of length 
dE/dz = — K, what is true both for the color string and for the bremsstrahlung [ pi 
mechanisms of hadronization. In this case we have 
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xi = 1 - '^Az , (22) 

where Az is the longitudinal distance covered by the cc pair in the color-octet state. 

The correction from this mechanism to the cross section for nuclear photoproduction 
can be calculated as follows 



7A^ °o 



A (^] = 27r ^ 5^ ^(7^^ ^ J/'^N) J d% J dz, pA{b,z^) J dz, p^i^z^) 



Z2- Zi-{1- Xi)- 



exp {-af/ im - T.,(6) + T.,(6)]} . (23) 



The derivation of this expression includes following ingredients. We assume that the photon 
energy is sufficiently high, u ^ M'^Ra/2 to allow the photon to convert into the cc pair 
long in advance. At the first interaction point {b,zi) the cc pair converts into the color 
octet state, corresponding to the cut Pomeron with cross section afj^. It is important 
that no longitudinal momentum is transferred to the cc pair at this point. All the loss of 
longitudinal momentum, (1 — Xi)z/, is due to the hadronization, until the last interaction 
at point {b,zi), where the color octet cc pair converts back to the colorless state. This 
interaction proceeds with the same cross section as the very first one, except for a factor 
1/8 (no summation over final colors of the cc in this case). We also use the relation cr*^ = 
{crf^fy/lGirBfi^, where Bfi^ ^ 5 GeV'"^ is the slope of \E'A^ elastic scattering (see section 
3. The 5-function in eq. (|23|) takes care of energy-conservation including the energy loss 
( |2^ ) during the intermediate state hadronization. We also assume that the cc pair does 
not attenuate when it propagates through the nucleus in a color octet state. Actually, 
rescatterings of such a pair also may cause an attenuation |2^, but it is small provided 



that the photon energy is sufficiently high (1 — Xi)v / n ^ Ra- However, the cc pair does 
attenuate propagating through the nucleus in the colorless state before the first interaction 
(if qeiRA ^ 1) and after the last one. We assume in eq. (|23|) a long coherence length and 
zero formation length. We neglect the energy variation due to these effects because the 
estimate eq. (p3|) is quite rough and the accuracy of data [|l], || is low as well. 
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The only unknown parameter we are left with in eq. (|23|) is k, the energy loss per unit 
length. In the string model it equals the string tension. The color-octet string tension k may 
substantially exceed the known value for color-triplet strings K3 ~ (27ra^)^^ ^ 1 GeV/fm 
| 23| , where a'j^ ~ 1 GeV^^ is the universal slope of Regge trajectories. Indeed, the string 
tension of the color-octet string must to be related to the slope of the Pomeron trajectory, 
a'p ~ 0.2 GeV^"^, which is much smaller than a^. Thus, ks ~ (27rap)^^ ^ 5 GeV/ fm 

One may also treat the energy loss perturbatively as a result of gluon bremsstrahlung by 
the color-octet cc pair. It was demonstrated in [^] that bremsstrahlung provides a constant 
density of energy loss, similar to the string model. The mean squared color-octet charge 
is 9/4 times bigger than the color-triplet one, so is the energy loss. Besides, in the case 
under discussion the color-charge appears during the first interaction with a bound nucleon 
at the time t = ti, then propagates for the time. At = t2 — ti and disappears when the 
cc pair becomes colorless. Such a process with switch-on and switch-off causes a double 
bremsstrahlung |^ (see also |^): only that part of the frequency spectrum is emitted 
during a time interval At, which has radiation time tr = 2uj/k^ < At, where u and k-r 
are the energy and transverse momentum of the radiated gluon, respectively. At the time 
t = t2 when the cc pair converts to a colorless state and the radiation process stops. A 
new radiation process caused by the charge stopping starts, and another piece of the gluon 
spectrum is radiated, identical to the one previously emitted. Thus, the intensity of radiation 
and the energy loss are twice as large as in a single-scattering process. Summarizing, we 
expect the energy loss for gluon bremsstrahlung in the double-scattering process with a 
color-octet intermediate state to be 4.5 times larger than that for radiation of a color-triplet 
charge, produced in a single scattering process. This estimate is very close to the one we 
found above in the string model, so we fix for further calculations k = 5 GeV/ fm. 

The new mechanism provides an xi-dependence of the cross section quite different from 
the standard one. There is no cross section for xi < 1 — kRa/ v-, since the amount of lost 
energy is restricted by the length of the path of the color-octet cc pair inside the nucleus. 
This contribution does not peak at Xi — > 1, but may have even a minimum due to a longer 
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path of the colorless cc pair in nuclear matter. 

The contribution of this mechanism may cause an A-dependence steeper than A^, since 
there is more room for integration over longitudinal coordinates Zi and Z2 ■ For this reason the 
contribution of the mechanism under discussion, is important for heavy nuclei at moderately 
high Xi- We compare the contribution of Glauber-Gribov mechanism calculated with eq. ( |T7| ) 
and the novel one given by eq. (p3D to J/"^ photoproduction on lead at different photon 
energies in Fig. 4. We see a strong energy variation of the latter contribution, which shifts 
with energy towards Xi = 1. At the same time the energy dependence of nuclear effects 
provided by a variation of the coherence and formation lengths is so small that we plotted 
only one curve for the Glauber-Gribov mechanism at z/ = 100 GeV. 

We see that both mechanisms the Glauber one and the two-step production are of the 
same order at Xi ~ 0.6 — 0.7. 

The novel mechanism proposed here does not display the Feynman scaling typical for the 
standard mechanism, but this contribution grows with energy at fixed Xi. This is easy to 
understand: the total contribution integrated over Xi is energy independent, but the range 
of 1 — Xi shrinks oc l/u. 

With the above value of k = 5 GeV/ fm and afj^ = 5 mb we calculate the contribution 
of the double-color-exchange mechanism (^). Combining with the Glauber contribution 
from section 2 we have 



The result for the ratio of the cross sections for tin to carbon is plotted in Fig. 3 for two 
photon energies u = 70 and 100 GeV typical for the NMC experiment in comparison with 
the data g. 

It should be noted that such a comparison is legitimate only for xi < 0.85 where the 
experimental resolution of the J/\l/ momentum suffices to separate events of inelastic pho- 
toproduction from the quasielastic ones. For xi > 0.9 quasielastic photoproduction of J/\l/ 
is supposed |1|, 0] to be the dominant process. In this region we should compare the data 




(24) 
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30 GeV 
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/ ' 
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0.4 0.6 0.8 1.0 
CK / 

Figure 4: The cross section per nucleon for J/\l/ photoproduction on lead. 
The sohd curve corresponds to Glauber-Gribov mechanism of nuclear 
shadowing at = 100 GeV . The novel mechanism contributions at 
V = 70, 100 and 140 GeV are shown by dashed curves. 



with the quasielastic photoproduction cross section. The corresponding formula is derived 
in 



I0| , |Tl|. It is different from formulas (0) - (|5|) for inelastic photoproduction, but in the 
approximation (j*f^(T) ^ 1 looks analogous to eq. ([T7|) 



ag^ii^A ^ vl>) ^ a(7iV ^ v^iV) A |l - ^afj (T) [l + [l + e R FHqj)] } , 

(25) 

The results of calculation with eq. (|T^ shown in Fig. 3 correspond to quasielastic scattering 
in the limit xi ^ 1, where the momentum transfer eq. (^) acquires the value qj = (M|, — 
M|)/2i/. 

Fig. 3 shows a good agreement with the data both for inelastic and quasielastic photo- 
production of J/\E' off nuclei. 
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It is important to study the energy dependence of the effect under discussion, which 
is expected to be interesting. It seems most hkely that new data on inelastic J/^ pho- 
toproduction will be taken with the HERMES spectrometer at HERA. In Fig. 5 we show 
predictions for nuclear suppression/enhancement factor S^"^ at u = 30 GeV for xenon and 
nitrogen. 
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Figure 5: Nuclear suppression/enhancement factor for J/^ photopro- 
duction on xenon and nitrogen at = 30 GeV 



7. Conclusions and discussion 

We summarize the main results of the paper as follows. 

• The cross section of inelastic photoproduction of charmonia off nuclei is calculated 
within the eikonal (Glauber) approximation. Corresponding formulas are derived for 
the first time. We also correct the Glauber formulas for inelastic shadowing, which 
makes nuclear medium substantially more transparent for charmonium. 

• The experimentally observed nuclear enhancement of inelastic J /'^ photoproduction 
at Xi < 0.85 is found to be too large compared with expectation of Glauber model to 
be explained by an enhancement of gluons in nuclei at low Xsj- 
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A novel QCD mechanism of charmonium photoproduction is suggested, which is not 
included in the standard multiple-scattering theory and can only occur in a nuclear 
target. This mechanism arises as a natural consequence of the colored structure of 
the Pomeron: the projectile hadron (or a hadronic fluctuation of the photon) induces 
a color dipole in the target nucleus, leaving the nucleus in a colorless state. This is a 
specific way of diffractive excitation of the nucleus and inclusive production of leading 
hadrons. 

As different from the Glauber model the new mechanism contribution is essentially 
energy-dependent and violates Feynman scaling. Due to steeper an A-dependence 
this mechanism successfully competes with Glauber one at moderate Xi in the energy 
range of the EMC and NMC experiments and nicely explains the observed [|l], |] 
nuclear enhancement of the photoproduction cross section. This effect is expected 
to be squeezed towards Xi = 1 and to vanish under the quasielastic peak at higher 
energies. 

We conclude that the observed nuclear enhancement of inelastic photoproduction of 
J/\l/ is the evidence of a new phenomenon related to specifics of QCD. An analogous 
effect in the case of interaction of light hadrons turns out to be suppressed and escapes 
observations [0, ^ . Photoproduction of heavy quarkonia might be a unique reaction 



where the new mechanism can be observed. 

Note, that a specific channel of decay of the color-dipole created inside the nucleus into 
two nucleons was suggested in as a source of the backward protons in the laboratory 



frame. Measurements by the SIGMA Collaboration of backward proton production in 
pion-beryllium interaction at 40 GeV nicely confirmed the theoretical expectations. This is 
another manifestation of the of double-color-exchange mechanism. 

The important signature of the mechanism under discussion its specific energy depen- 
dence. The nuclear enhancement does not vanish at higher energies, but is squeezed towards 
Xi = 1. This is different from gluon distribution in a nucleus, which anyway have to switch 
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to nuclear shadowing at higher energies, i.e. smaller Bjorken X2- 
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Appendix 

We estimate the correction S^dal"^ /dxi) from inelastic production of J/\E' with inelastic- 
ity a at one point in the nucleus, with subsequent inelastic rescattering with inelasticity (3 
on another bound nucleon. This contribution to the cross section integrated over transverse 
momenta reads 

S ( ^] ^Af Ida jd/3 - a/?) ^ ^ , (26) 
\ dx-i I J J da dp 

where 

T = ^Jd'b T\h) exp[-af/T(6)] (27) 
We assume for the sake of simplicity that the energy is high, qeiRA ^ 1- 



After integration (26) takes the form 



V dxi I dxi 



noi where the factor K{xi) is a slow function of xi. 



K(xi) 



'afr\ Gjpp ^ pX) 



[1 + Xi) In ( ^-^ + Xi In ^ — 



(29) 



Here e = M^/s, where Mq is the minimal mass of the excited nucleon in final state, which 
we include in the integration over inelasticities in (|26|). At Mq = 2 GeV and xi = 0.9 we 
estimate the factor K ^ 0.03 fm^ 
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The effective nuclear thickness (27) for tin is T'^" ~ 0.3 fm~'^. 

We conclude that according to these estimates and expression (|29D the correction under 
discussion has nearly the same form of Xi-dependence as the standard mechanism contri- 
bution, but is suppressed by the factor K T A/S^^ ~ 0.016. We can neglect such a small 
correction for our purposes. 
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